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The ground and excited state properties of small helium clusters, ^Hcn, con- 
taining nanoscale (^3-10 A) planar aromatic molecules have been studied 
with quantum Monte Carlo methods. Ground state structures and energies 
are obtained from importance- sampled, rigid-body diffusion Monte Carlo. 
Excited state energies due to helium vibrational motion are evaluated us- 
ing the projection operator, imaginary time spectral evolution technique. We 
examine the adsorption of N helium atoms (N < 24) on a series of planar 
aromatic molecules (benzene, naphthalene, anthracene, tetracene, phthalo- 
cyanine). The first layer of helium atoms is well-localized on the molecule 
surface, and we find well-defined localized excitations due to in-plane vibra- 
tional motion of helium on the molecule surface. We discuss the implications 
of these confined excitations for the molecule spectroscopy. 

PACS numbers: 36.40.Mr, 6l.46.-hw, 67.40. Yv, 67.70. -hn. 

1. INTRODUCTION 

The utiUty of He II as a gentle, quantum matrix for the high-resolution 
spectroscopic study of impurity molecules is now well-established.— Accom- 
panying these recent developments in experimental methodology is an in- 
terest in examining the spectroscopy of large, planar aromatic molecules 
(PAM) in helium. These PAMs can be viewed as nanoscale precursors to 
bulk graphite surfaces, whose size and geometry can be systematically tuned. 
The experiments have examined the electronic spectra of the molecules, and 
can be divided into two general classes: a) those involving small numbers of 
helium atoms {N < 17) around the impurity molecule, and b) those where 
the PAM is embedded in a large helium droplet (A^ ~ 10000). 
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The electronic spectra of PAMs in helium exhibit features not found for 
the corresponding bare molecules. These new features are presumably due 
to the presence of helium on the molecule surface, but a more specific under- 
standing of their origins does not yet exist. In the absorption experiments, 
sharp excitations have been observed at ~ 3 — 30 K above the electronic ori- 
gin, which have been qualitatively attributed to the localized vibrations of 
helium atoms adsorbed on the molecule surface. ^'^'^ This energy range also 
encompasses the band of states (phonon wings) associated with collective 
compressional excitations of the cluster interior,>5. and in the large cluster 
experiments these localized excitations appear as additional structure su- 
perimposed on the phonon wing sidebands. 

For the electronic spectra of a number of different PAMs (tetracene, 
Mg-phthalocyanine, and various indole derivatives) in large He droplets, the 
peak associated with the electronic origin (zero-phonon line) also exhibits 
additional structure, with splittings on the order of ~ 1.5 K.-Sii These split- 
tings are not observed in the small cluster experiments, - and their physical 
origins are not completely understood at this point. One possibility is that 
they are due to the excitation of low-energy modes involving the collective 
motion of helium atoms on the molecule surface. Theoretical calculations for 
the HeAT-benzene cluster suggest that individual helium atoms are strongly 
localized on the molecule surfacc^^^i As we discuss further below, this is a 
general feature which we find for a number of other PAMs, so another pos- 
sible explanation is that these splittings do not arise from elastic transitions 
at all, but are instead due to inhomogeneities in the local helium environ- 
ment around the molecule. Recent emission experiments have also revealed 
additional features not previously seen in absorption. For a number of other 
PAMs (phthalocyanine, Mg-phthalocyanine, pentacene), when an electronic 
and vibrational mode in the impurity molecule are simultaneously excited, 
the resulting emission spectrum as the molecule relaxes to the ground state 
is similar to that for absorption, except each line is split into two peaks that 
are separated by ~ 14.8 K.-^ This suggests that the helium dynamics after 
electronic excitation are more complicated than previously thought. 

In this work, we present the results of quantum Monte Carlo calculations 
for small numbers of helium atoms {N < 24) around various PAMs. The 
strength and anisotropy of the molecule-helium interaction gives rise to a 
strongly localized layer of helium atoms on the molecule surface in the cluster 
ground states. This localization has important consequences for the excited 
states,— and we report some preliminary results for excitations due to the 
vibrational motion of helium atoms in this first monolayer. We propose 
that the strong ground state localization implies similar confinement of the 
excited state, which may consequently affect the dynamics of energy transfer. 



structure and Energetics of Helium Adsorption on Nanosurfaces 



2. THEORY AND METHODS 

The Hamiltonian for the ^Hcat-PAM cluster consists of the rigid-body 
kinetic energy for the PAM, the translational kinetic energy for the ^He 
atoms, and the total potential energy V. This V is taken as a sum over 
two-body helium-helium— and helium-PAM interactions, where the helium- 
PAM interaction V[ is a sum of atom-atom Lennard-Jones pair potentials in 
the principal axis frame of the molecule: 



yi(r)=^4e. 



12 



(1) 



The sum over a = C,H,N runs over the individual atoms of the molecule 
situated at the position r^, and Lennard-Jones parameters e^, cTq were taken 
from Vidali et al. The potentials evaluated in this manner are qualitatively 
similar. Each of the potentials for the linear PAMs exhibits minima near 
the center of the aromatic rings. The global minima in the potential of 
each molecule occur at ~ 2.7 A above and below the molecule surface. The 
potential for tetracene at this distance from the molecule surface is shown 
in Figure 

Ground state properties are obtained using diffusion Monte Carlo (DMC) 
methods, in particular the importance-sampled rigid body diffusion Monte 
Carlo algorithm.^ The trial functions used here are products of two-body 
factors, 

N N 

= n n e-*H-(i''^--'''i) (2) 

j=l i<j 

where ti and tpje describe helium-PAM and helium-helium correlations, re- 
spectively: 

t,(r) = Y,(r^y +^^' + by' + cz', tMr) = (—y. (3) 

The trial function parameters Cq,, CHej b, c are obtained by minimizing the 
variational energy of H with respect to ^t- Excited state energies are 
obtained using the projection operator, imaginary time spectral evolution 
(POITSE) method.— This yields the spectral density function 

N 

KiE)cxJ2\i^n\A^\-^T)\^^{E-Er. + Eo), i^ = E«]' W 

n j=l 

where {(/>n} and {En} are the set of eigenfunctions and associated eigenval- 
ues of the Hamiltonian H, respectively. The many-body operator is a 
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Fig. 1. Helium density distributions at a distance of z = 2.9 A above the 
molecule plane for (a) ^He4-naphtalene, (b) ^Heg-anthracene, and (d) ^Heg- 
tetracene. A cut of the tetracene potential at z = 2.7 A is also shown in (c), 
where contour lines run from —160 K to K, in increments in 15 K. The 
tetracene molecule on the z = plane is drawn in bold lines. 

Bose-symmetrized sum of one-body excitation operators , and is chosen to 
connect the ground state to the excited state(s) of interest. 

3. RESULTS AND DISCUSSION 

Ground state helium density distributions for ^He4-naphthalene, ^Hee- 
anthracene, ^Heg-tetracene, and ^He24-phthalocyanine were obtained using 
the quantum Monte Carlo methods described above. The number of ^He 
atoms around each of the PAMs studied here corresponds to full coverage 
of both sides of the nanosubstrate surface. In all cases, the global density 
maxima occur at a distance of z ~ 2.9 — 3.2 A above and below the molecule 
plane. Figures ^p, and show cuts of the local helium density at 
this distance, along a plane parallel to the molecule surface. It is apparent 
that the local helium densities consist of well-separated peaks, each corre- 
sponding to a single ^He atom. For the linear PAMs, each individual ^He 
atom is approximately situated above a six-membered carbon ring. For ph- 
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Fig. 2. Left: Helium density distribution at a distance of z = 3.2 A above 
the molecule plane for ^He24-phthalocyanine. The phthalocyanine molecule 
on the z = plane is drawn in bold lines. Right: Spectral density function 
for ^Hes-tetracene. 



thalocyanine, we find that this surface can support up to 12 He atoms on 
each side. At full coverage, this 2D layer of '^He on phthalocyanine consists 
of four "^He atoms grouped near the center of the molecule, surrounded by 
an outer ring of eight ^He atoms. 

Previous work with the ^HeAr-benzene system (iV = 1,2,3, 14) has re- 
vealed similar ground state features.^ The benzene system can be viewed as 
the simplest ^Hctv-PAM system, serving as a building block for models of 
more complicated one- and two-dimensional nanosubstrates. In that study, 
not only were the ground states observed to be strongly localized, but the ex- 
cited states also exhibited localized character. That is, excitations due to the 
localized vibrational motion of ^He on the benzene surface were found, with 
energies of up to ~ 23 K above the ground state. We propose here that these 
excited state features exist generally for larger '^He-PAM systems, and that 
they are a result of the strong localization of the helium ground state den- 
sity distribution. Figure shows the spectral density function [Eq. (jlj)] for 
^Hes-tetracene, evaluated using the POITSE methodology discussed above. 
For a choice of = x and x^, where the rc-direction is taken to be along the 
long axis of the molecule, we find excitations at 6.3, 9.9, 20.8 K above the 
ground state. These states correspond to the collective vibrational motion of 
the ^He atoms moving on the molecule surface along the x-direction. Work 
is currently in progress to extend these calculations to larger A^, in order 
to ascertain whether the collective character of these excitations persists as 
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additional atoms are added around the molecule. 

A helium excitation confined on the molecule surface would have impor- 
tant consequences for spectroscopy, particularly for the dynamics of energy 
dissipation. Typically, the molecule is excited from the electronic ground 
state Sq to an electronic excited state Si. For a weakly interacting small 
molecule in a vibrationally excited state of Si , this excess vibrational energy 
is transferred to the local helium environment. It is generally assumed that 
this energy is rapidly carried away from the molecule due to the large ther- 
mal conductivity of He II.— In a helium droplet, this energy is dissipated 
as helium atoms boil off the surface of the droplet, and so the relaxation 
back to Sq derives primarily from the ground vibrational state of . On the 
other hand, if the excess vibrational energy of the molecule is transferred 
to the helium and trapped near the surface of the molecule, the resulting 
local helium environment could be very different from that of the ground 
vibrational state and would then correspond to a metastable state. These 
two situations could give rise to two distinct Si surfaces, resulting in a split 
emission spectrum similar to that observed by Lehnig and Slenczka.^ Calcu- 
lation of the local helium density for the Si state, which requires a currently 
unknown interaction potential, would facilitate a more detailed analysis. 
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